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Abstract 

In this second of three articles on the optical absorption of electrons in a half- 
filled Peierls-distorted chain we present exact results for strongly correlated 
tight-binding electrons. In the limit of a strong on-site interaction U we map 
the Hubbard model onto the Harris-Lange model which can be solved exactly 
in one dimension in terms of spinless fermions for the charge excitations. The 
exact solution allows for an interpretation of the charge dynamics in terms 
of parallel Hubbard bands with a free-electron dispersion of band-width W, 
separated by the Hubbard interaction U. The spin degrees of freedom enter 
the expressions for the optical absorption only via a momentum dependent 
but static ground state expectation value. The remaining spin problem can be 
traced out exactly since the eigenstates of the Harris-Lange model are spin- 
degenerate. This corresponds to the Hubbard model at temperatures large 
compared to the spin exchange energy. 
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Explicit results are given for the optical absorption in the presence of a 
lattice distortion 5 and a nearest-neighbor interaction V. We find that the 
optical absorption for V = is dominated by a peak at u = U and broad but 
weak absorption bands for \u — U\ < W. For an appreciable nearest-neighbor 
interaction, V > W/2, almost all spectral weight is transferred to Simpson's 
exciton band which is eventually Peierls-split. 
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I. INTRODUCTION 



Some charge-transfer salts are understood as strongly correlated one- dimensional elec- 
tron systems for half-filled bands (Mott-Hubbard insulators) (Farges 1994); (Alcacer, Brau, 
and Farges 1994). The extended Hubbard model for interacting electrons on a Peierls- 
distorted chain at half-filling is considered appropriate for these materials (Mazumdar and 
Dixit 1986); (Fritsch and Ducasse 1991); (Mila 1995). There are only few studies of the 
optical, i. e., finite frequency properties of correlated electron systems since their calculation 
is a formidable task (Kohn 1964); (Maldague 1977); (Lyo and Galinar 1977); (Lyo 1978); 
(Galinar 1979); (Campbell, Gammel, and Loh 1988); (Mahan 1990); (Shastry and Suther- 
land 1990); (Stafford, Millis, and Shastry 1991); (Fye, Martins, Scalapino, Wagner, and 
Hanke 1992); (Stafford and Millis 1993). 

In this second article on the optical absorption of electrons in half-filled Peierls-distorted 
chains we present a detailed analysis of the optical absorption in the limit of strong corre- 
lations and for a half-filled band where the charge and spin dynamics decouple. We find 
that the physics of the half-filled Hubbard model at strong correlations is determined by 
the upper and lower Hubbard band for the charges which are parallel bands. This essential 
feature and its significant consequences have been missed in earlier analytical and numerical 
investigations. In this work we include a finite lattice dimerization and nearest-neighbor 
interaction between the electrons, i. e., we analyze the extended dimerized Hubbard model 
at strong correlations. 

The paper is organized as follows. In section [H] we address the Hubbard model at strong 
coupling from which we derive the Harris-Lange model which determines the motion of the 
charge degrees of freedom. In the ground state there are no free charges but only singly 
occupied lattice sites. The exact spectrum and eigenstates of the Harris-Lange model are 
presented in section |H| for the translational invariant and the dimerized case. The results can 
be interpreted in terms of two parallel Hubbard bands for the charges which are eventually 
split into Peierls subbands. Optical absorption can now formally be treated as if we had 
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independent (spinless) Fermions. 

Unfortunately, this simple band structure interpretation is blurred by the spin degrees of 
freedom which enter the expressions for the optical absorption in terms of a very complicated 
ground state expectation value. In section [TV] we treat the case of the Harris-Lange model 
where all spin configurations are equally possible ground states. This corresponds to the 
Hubbard model at temperatures large compared to the spin exchange energy. It allows 
the calculation of the optical absorption even in the presence of a Peierls distortion and 
a nearest-neighbor interaction between the charges. A summary and outlook closes our 
presentation. Some details of the calculations are left to the appendices. 



II. STRONGLY CORRELATED MOTT-HUBBARD INSULATORS 

A. Tight-binding electrons on a Peierls-distorted chain 

For narrow-band materials the electron transfer is limited to nearest neighbors only. In 
standard notation of second quantization the Hamiltonian for electrons in the tight-binding 
approximation reads 

T(S) = -t jr (l + (-l)'tf) (c+ c l+ha + c+ (1) 

l=l,a 

where 5 describes the effect of bond-length alternation on the electron transfer amplitudes. 
As usual the Hamiltonian can be diagonalized in momentum space. We apply periodic 
boundary conditions, and introduce the Fourier transformed electron operators as c\ a = 
\A/-k Td=i exp(ikla)cl >a for the L momenta k = 2nm/(La), m = -(L/2), . . . (L/2) - 1. We 
may thus write 

|/fc|<7r/(2a),cr 

with the dispersion relation e(k) and hybridization function A(k) defined as 

e(k) = —2tcos(ka) (3a) 
A(k) = 2t5 sm(ka) . (3b) 
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The Hamiltonian can easily be diagonalized in /c-space. The result is (Gebhard, Bott, 
Scheidler, Thomas, and Koch I 1996) 



f (5) = £ E(k)( 



a k,a,+ a k,a,+ a k,a,- a k,a,- 



(4) 



\k\<Tr/(2a),cr 

Here, ±E(k) is the dispersion relation for the upper (+) and lower (— ) Peierls band, 



E(k) = Je(k) 2 + A(k) 2 . 



(5) 



The new Fermion quasi-particle operators ajj" ± for these two bands are related to the 
original electron operators by 



a Kc _ = a k c k(7 + i(3 k c k +ir,* 

^k,a,+ = PkC kt<7 — ia k C k+7V(J 



(6a) 
(6b) 



with 



ot-k 



\2\ E{k)j 



(3 k = 



1 



1 



e(k) 



|sgn(A(fc)) 



(7a) 



(7b) 



\ 2 V 1 E ( k ), 

Details of the upper transformation and the optical absorption of this model are presented 
in (Gebhard et al. I 1996). 



B. Hubbard model 

The only spinful interacting electron model that can be solved exactly for all values of the 
interaction strength is the Hubbard model in one dimension (Hubbard 1963); (Gebhard and 
Ruckenstein 1992); (EBler and Korepin 1994); (Gebhard, Girndt, and Ruckenstein 1994); 
(Bares and Gebhard 1995). For narrow-band materials the electron transfer is limited to 
nearest neighbors only, and the interaction is supposed to be described by the purely local 
(Hubbard-) interaction of strength U, 
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-^Hubbard — T + UD 

^E^Evu, (8) 
i i 

where ri; j(J = cf a c la . is the local density of cx-electrons, and T = T{5 = 0). 

The model (H) poses a very difficult many-body problem. Its spectrum and, in particular, 
its elementary excitations can be obtained from the Bethe Ansatz solution (Lieb and Wu 
1968); (Shastry, Jha, and Singh 1985); (Andrei 1995). Its low-energy properties including the 
DC-conductivity, <tdc = Re{a(u; = 0)}, can explicitly be obtained from the corresponding 
g-ology Hamiltonian (Schulz 1990); (Schulz 1991) or from conformal field theory (Frahm and 
Korepin 1990); (Frahm and Korepin 1991); (Kawakami and Yang 1990); (Kawakami and 
Yang 1991). The Hubbard model describes a (correlated) metal for all U > for less than 
half-filling. Unfortunately, the Bethe Ansatz solution does not allow the direct calculation 
of transport properties at finite frequencies. 

At half-filling the one-dimensional Hubbard model describes a Mott-insulator which im- 
plies that o"dc — for all U > 0. The density of states for charge excitations displays two 
bands, the upper and lower Hubbard band, separated by the Mott-Hubbard gap. This gap 
is defined as the jump in the chemical potential at half filling, 

A MH = /i + (iV = L)-/i-(iV = L) 

= [E(N = L + 1) - E(N = L)] - [E(N = L) - E(N = L - 1)] . (9a) 

As shown by Ovchinnicov (Ovchinnicov 1969) the Mott-Hubbard gap can be obtained from 
the Lieb-Wu solution (Lieb et al. 1968) in the form 



A MH = -T7- \ . , - /m 9b 
U Ji smh{27Tty / U J 



(2W/ir)J4U/Wexp{-irW/{2U)) for U < W = At 

(9c) 

U-W+ \n(2)W 2 /(2U) + 0(W 3 /U 2 ) for U > W = U 
It is obvious that optical absorption is only possible if u > A MH . It is further seen that the 
upper and lower Hubbard band are well separated for U 3> W . 
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One might expect that the optical absorption for large interactions, U ^> W, and high 
temperatures, k B T 3> J = 0(W 2 /U), shows the signature of a broad band-to-band transi- 
tion for U — W < u < U + W (units H = 1), similar to the Peierls insulator (Gebhard et 
al. I 1996). Such considerations seemed to be supported by analytical (Lyo et al. 1977); 
(Lyo 1978); (Galinar 1979) and numerical calculations (Campbell, Gammel, and Loh 1989). 
Below we will calculate a{uj > 0) in the limit U 3> W, and show that the linear absorption 
is actually dominated by a singular contribution at uj = U because the upper and lower 
Hubbard band are in fact parallel bands. The situation changes for k B T <C J which we will 
consider in (Gebhard, Bott, Scheidler, Thomas, and Koch III 1996). 

C. Harris-Lange model 

In the following we will address the limit W/U — > where matters considerably simplify 
since the charge and spin degrees of freedom completely decouple (Ogata and Shiba 1990); 
(Parola and Sorella 1990). For example, for less than half- filling, N < L, and U — oo the 
eigenenergies become those of a Fermi gas of N h — L — N holes with dispersion e(fc), and 
each energy level is 2 iV -fold degenerate in the thermodynamical limit (Beni, Pincus, and 
Holstein 1973); (Klein 1973); (Ogata et al. 1990); (Parola et al. 1990). 

To facilitate the discussion of the strong coupling limit we map the Hubbard model 
onto a problem for which the number of double occupancies is conserved. For a large on-site 
Coulomb repulsion W/U — > it is natural to start with a spectral decomposition of operators 
into those which solely act in the upper or lower Hubbard band, and to perturbatively 
eliminate those parts in -^Hubbard which couple the two bands. For the Hubbard model 
this has first been achieved by Harris and Lange (Harris and Lange 1967); (van Dongen 
1994), and the resulting effective Hamiltonian to lowest order in W/U will thus be called 
the "Harris-Lange" model. It offers several advantages, both for analytical and numerical 
calculations. 

To carry out the spectral decomposition we start from the case t — 0. The Fermi 
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annihilation operator can be split into a part which destroys an electron on a single occupied 
site and does not change the energy of the state, and another part which destroys an electron 
on a double occupied site and thus decreases the energy by U, 

ci, <T = n l - IT c l><T + (l-n l - <T )c lta . (10) 

The corresponding creation operator can be treated accordingly. 

If we now turn on the hopping of electrons (t 7^ 0) we may split the kinetic energy 
operator into 

f = f LHB + T UH B + T+ + f (11a) 

Tlhb = (-*) £(i- hi,-„) (cTA+i^ + J C 1 - ( llb ) 

1,(7 

TUHB = H) £ hx-a (ctA+l,a + 4+1^) ^1+1 -a (He) 
l,a 

T + = (~t) £ [ni,-act,A+i, ff (1 - n i+1 - a ) + n i+ i,_ CT c z + l (7 c l a (1 - nj (lid) 

l,a 

f- = (f+Y . (lie) 

The operator Tlhb for the lower Hubbard band describes the hopping of holes while doubly 
occupied sites can move in the upper Hubbard band via Tuhb- Their number is conserved 
by both hopping processes. These two bands will constitute the basis for our approach. The 
operator T + (T~) increases (decreases) the number of double occupancies by one. 

Similar to the Foldy-Wouthuysen transformation for the Dirac equation (Bjorken and 
Drell 1964) we apply a canonical transformation that eliminates the operators T^ 1 to a given 
order in t/U, 

%o = e^c^e-^ (12) 

with (S'(c)) = S(c) . As shown by Harris and Lange (Harris et al. 1967); (van Dongen 
1994) the operator to lowest order in t/U reads 



which can easily be verified since 



D,T 



± 



= ±T±. 

The transformed Hamilton operator in the new Fermions becomes the Harris-Lange 
model 

HL = f[ HB + T£ HB + UD C , (14) 

if we neglect all correction terms to order t/U and higher. The energies obtained from 
the Harris-Lange model thus agree with those of the Hubbard model to order t(t/U)~ l 
and t(t/U)°. For all other physical operators which do not contain a factor of U/t we may 
replace 

Cl,a = . (15) 

because the error is only of order (t/U). In the following we will thus make no distinction 
between the operators q j(J and q )CT to lowest order in t/U. 

The Hamiltonian has the following symmetry. The particle-hole transformation 

%^t T ^ = iKe M \_ a (16a) 

"W^ph 1 = %\cKla-k,-a ( 16b ) 

with Aj = — = 1 is generated with the help of 

r ph = e™i^++c-) = n [i - ( A + m + + en] (iec) 

C + =pj + = EQ + = D-l)% ; ^ = (l-n u )(l-n u ). (16d) 
I I 

The additional phase factors i\ a are irrelevant global phases, and can be ignored since 
there is always an equal number of Fermion creation and annihilation operators of each spin 
species. The operators for the motion of holes and double occupancies are mapped into each 
other, 

7\jHB | — ► TlhB TLhB ^ TuHB • (17) 
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Furthermore, T UH b + Tlhb, C <± = 0. This symmetry allows for an exact solution of the 
model since there is essentially no difference in the motion of double occupancies in the 
upper Hubbard band and holes in the lower Hubbard band. 

The discussion above is readily generalized to the case of dimerization in the Harris-Lange 
model. The model Hamiltonian reads 

Hhl — Tlhb(S) + T UH b (5) + UD (18) 
in an obvious generalization of the kinetic operators for the upper and lower Hubbard bands. 

D. Optical absorption and optical conductivity 

The dielectric function e{uj) and the coefficient for the linear optical absorption a(u) are 
given by (Haug and Koch 1990) 

e » = 1 + (19a) 
to 

n b c 

where Re{. . .} denotes the real part and n& is the background refractive index. It is sup- 
posed to be frequency independent near a resonance. Hence, the real part of the optical 
conductivity directly gives the absorption spectrum of the system. 

The standard result (Maldague 1977); (Mahan 1990) for the real part of the optical 
conductivity in terms of the current-current correlation function x{°°) is 

Re{a(.)} = (20 ) 
to 

xM= lT/o ^MtWJ (21) 

where Af± is the number of chains per unit area perpendicular to the chain direction. 

The current-current correlation function can be spectrally decomposed in terms of exact 
eigenstates of the system as 
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± 



La 



EI(o|j»I 2 



i 



i 



(22) 



u + (E n - E Q ) + 17 w - (J5„ - £ ) + «7_ 
Here, |0) is the exact ground state (energy i? ), \ n ) are exact excited states (energy E n ), and 
|(0|j|n)| 2 are the oscillator strengths for optical transitions between them. Although 7 = + 
is infinitesimal we may introduce 7 > as a phenomenological broadening of the resonances 
at oj = ±(E n — Eq). The spectral decomposition of the real part of the optical conductivity 
reads 

Af\_n 



Re{a(uj)} 



EK°b»l P (w - (E n - E Q )) -5(u+ (E n - E ))] 



Lauj ,, 

which is positive for all u. 

In the following we will always plot the dimensionless reduced optical conductivity 

uRelaiuj > 0)} 

a Ted {u > = 1 1 — ,s . 

M±ae 2 W 

Furthermore we replace the energy conservation 5(x) by the smeared function 



(23) 



(24) 



6(x) 
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(25) 



ir(x 2 + 7 2 ) 

to include effects of phonons, and experimental resolution. For all figures we graphically 
checked that the sum rules of appendix |A 1| were fulfilled. 



E. Current operator 

As derived in (Gebhard et al. I 1996) the current operator is given by 

j=-eJ2 if a (l + (l + (-l)S) (St+iA* ~ ^ g J+v) ( 26 ) 

La 

where rj = —\Ri + i — Ri — a\/a < is the relative change of lattice distances due to the Peierls 
distortion. Note that S and rj always have opposite sign. The current operator can be split 
into two parts, j = j^ tra + i^ter; where j^ tra moves electrons between neighboring sites 
without changing the number of double occupancies or holes. This (Hubbard-) intraband 
current does not change the number of double occupancies. Hence it can be ignored for the 
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optical absorption in the Harris-Lange model at half-filling. The current operator between 
the two Hubbard bands j^ tcr can be written as 



Winter Winter, + Winter, 



(27a) 



jL*, + = "(^) £ (l + (-1)'*) (f + (-1)'„) 

l,a 

-{itea) £ (l + (-1)'*) + 

l,a 



(27b) 



7 H 

J inter,— 



(27c) 



[(1 - h l+1 _ CT ) c^. 1)tr c J)(T ni _ ff - (f - nj c^q +li(T n i+ i _ 

where J^ t er,± create and destroy a neighboring pair of double occupancy and hole, respec- 
tively. 

Next we study the action of j^ tcr + on a pair of neighboring spins in a state in position 
space. It is a sequence of singly occupied sites (<r), holes (o), and double occupancies (•) 
from site 1 to L, e.g., 



|*> = |Tl,«2, 


°3,°4, J.5, 


• • Tl-3) *L-2, *l) 


(28) 


We introduce the notations 








• • • , (Tz, ± 1;, Ti+i) 


,...) = |. 


• , tj) • • •) =t • • • , Ui ii+i, ■ ■ ■) 


(29a) 


= 1, Sf jl+1 


= 1> = I- 


■ , ii, ii+i, ■ ■ ■) 


(29b) 


= 1, »Sm+i 


= 0) = |. 


■ , (ii, ii+i + u, ii+i), ■ ■ •) 


(29c) 


S«,«+i = 1, Sij+i = 


-1> = I- 


• , 1/) u+i, ■ ■ ■) 


(29d) 


= o, s* l+1 


= 0) = |. 


■ , (ii, ii+i - ii, ii+i), ■ ■ •) 


(29e) 



as the local spin triplet and spin singlet states. Furthermore, 

|C M+ i = 1, C* l+1 = 0) = | (•,, o l+1 - o h m l+1 ), . . .) 



(29f) 



denotes the local charge triplet state since C + \Cij+i = l,C^ l+1 = 0) ^ 0. With these 
definitions one finds 
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iinter,+ |S M+ i = 1) = (30a) 

j£ter,+ |S M+ i = 0) = -itea{\ + (-1)^)(1 + (-l) Z 7 7 )(-2)|C M+ i = 1, q m = 0) . (30b) 

It is thus seen that j^ ter preserves the spin of a neighboring pair such that AS = AS Z = 
is the selection rule for the spin sector. The selection rule for the charge sector is found as 
AC = 1, AC Z = 0. Note that the current operator does not commute with C as defined 
in eq. ( |16d|) . 

Finally, the current operator is invariant against translations by one unit cell and thus 
preserves the total momentum modulo a reciprocal lattice vector (Q = 1-nja for 5 = 0, 
Q = it /a for 5^0). However, the current operator can create or destroy a charge excitation 
with momentum q, and create or destroy a spin excitation with momentum — q. Although 
there is charge-spin separation in the Hubbard model for strong coupling the current operator 
mixes both degrees of freedom. This renders the calculation of the optical absorption of the 
Hubbard model a very difficult problem even in the limit of strong correlations. 

III. EXACT SOLUTION OF THE HARRIS-LANGE MODEL 
A. Translational invariant case 

The Harris-Lange model can exactly be solved by an explicit construction of all eigen- 
states. This has recently been shown by (de Boer, Korepin, and Schadschneider 1995) and 
(Schadschneider 1995) for periodic, and by (Aligia and Arrachea 1994) for open bound- 
ary conditions. Since optical excitations conserve total momentum we work with periodic 
boundary conditions where the total momentum is a good quantum number. 

The number N$ of sites with spin (singly occupied sites), and the number Nq = N^ + Nh 
of sites with charge (double occupancies and holes) are separately conserved in the Harris- 
Lange model. We have Nc + Ns = L lattice sites, N = Ns + 2Na electrons, and choose L to 
be even such that our lattice is bipartite for all L. In the sequence of singly occupied sites, 
double occupancies and holes of the state in eq. (|28|) we may identify subsequences for 
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the spins and the charges only (independent of the position on a special site). Additionally, 
the indices lj indicate the actual position of the charges Cj. The positions occupied by the 
spins are then the ones left over by the charges: 

1^) = \(h, h, ■ ■ ■ Inc-u Inc)'i (Ci, C 2 , • • • CjVc-i) Cn c )', (Si, S 2 , ■ ■ ■ Sjv s -ij Sn s )) (31a) 
where in our example 

(S 1 ,S 2 ,...S Ng - U S Ns ) = (tU--- J [A) (31b) 
is the subsequence for the spins and 

(Ci, C 2 , C 3 , . . . Cjv _i, CW ) = (•, o, o, . . . •, •) (31c) 
for the charges. The sequence for the positions occupied by charges is 

(h, l 2 , l 3 ,...l Nc -l, l Nc ) = (2, 3, 4, ... L — 2, L) . (31d) 

Since there is nearest-neighbor hopping only both the spin and charge subsequences are 
separately conserved up to cyclic permutations due to the periodic boundary conditions. 

To include the boundary effect we follow (de Boer et al. 1995) and (Schadschneider 1995) 
and introduce the properly symmetrized spin and charge sequences. To this end we define 
the operator for a cyclic permutation of the spin sequence, 

1s(Si,S2,-.-Sn s ) = (Sn s ,Si,...Sn s -i) , (32) 

and, equivalently, T c f° r the charge sequence. Let K s and K c be the smallest positive 
integers such that {t~s) S an d (^c) ° act as identity operators on a given spin and charge 
sequence, respectively. Then we define the K S K C states 

\Jk s K c \(1\, • • • In c )] (Ci, . . . C Nc ) kc ; (Si, . . . S Ns ) ks ) = (33) 

e E^A K jr K jr\* v s*s^W7Z*7Z°\(! 1 , . . . l Nc ); (Ci, . . . C Nc ); (Si, . . . S Ns )) 

^s=0 1^0=0 
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with the momentum shifts ks = litmsl (Kg a), ms = 0, 1, . . . (Ks — 1), kc = l-amcj (Kca), 
mc = 0, 1, . . . (Kc — 1). An extra phase factor (— for each double occupancy at site l 3 - has 
been included here through the operator exp(J]z vkIDi). This allows to make direct contact 
to the model considered by (de Boer et al. 1995); (Schadschneider 1995), and (Aligia et 
al. 1994). There the hopping amplitudes for lower and upper Hubbard band have opposite 
signs. 

We transform the states of eq. (|33|) into momentum space with respect to the charge 
coordinates. The exact eigenstates can now be classified according to their spin and charge 
sequence, their momentum shifts kc-, ks, and Nc momenta from the set of kj = 2nmj/ (La), 
rrij = — (L/2), . . . (L/2) — 1. The normalized eigenstates read 

L N °/ 2 \k h ...k Nc , (C h . . . C Nc ) kc ; (S 1 , . . . S Ns ) ks ) = (34) 
Hi- 1 ) 7 * ex P ( ia l r(j) ( k J + $ cs)) • • • In c ); (Ci, • • • C Nc ) kc ; (S u . . . S Ns ) ks ) 

h<...<l Nc V \ j=l J 

where the permutations V generate a simple Slater determinant for the momenta and the 
positions of the Nc charges, and $cs = (kc — ks) / L is an additional momentum shift which 
vanishes in the thermodynamical limit. 

It is straightforward (de Boer et al. 1995); (Schadschneider 1995) but lengthy to explicitly 
show that the states in eq. fl34|) have energy and momentum 

N c 

E = Y J <k 3 + $cs)+UN d (35a) 

3=1 

N C 

P = k s + (n /a)(N d -l) + Y,( k J + ®cs) + (N mod 2)tt/ a mod 2n/a (35b) 

3=1 



with e(k) given by equation (paf ). The essential arguments are repeated in the appendices |A 2 
and 

The above set of eigenstates is complete. After summing over the subspaces with differ- 
ent Ks, Kc the number of states which become degenerate in energy in the thermodynam- 
ical limit is given by 2 Ns Ncl/(NJNh\). The number of possible choices for the momenta 
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is L\/(Nc\{L — Nc)l) since the momenta are those of a gas of spinless Fermions. Altogether 
one finds for an even number of electrons N 

where eq. (10.33.5) of (Hansen 1975) and eq. (22.2.3) of (Abramovitz and Stegun 1970) have 
been used. This exhausts the Hilbert space for fixed number of electrons N. 



B. Dimerized Harris-Lange model 

The Harris-Lange model can also be solved in the presence of a finite lattice distortion 
as long as there is hopping between nearest neighbors only. For the U = oo Hubbard model 
at less than half filling this has already been realized some time ago (Bernasconi, Rice, 
Schneider, and StraBler 1975). The exact eigenenergies are those of spinless Fermions on a 
dimerized chain. 

For the Harris-Lange model with Nc charge excitations we choose momenta \ kj\ < n /(2a) 
of the reduced Brillouin zone, and one of the 2 Nc sequences (r 1; . . . t Nc ) with Tj = ±1. Let 
us introduce the operator IT r by 

U r \kx, . . .k r , . . . fejv c ; (Ci, . . . CN c )k c ', (Si, . . . SN s )k s ) 

= \ki, ...kr + Tv/a,... k Nc ; (d, . . . C Nc ) kc ; (Si, • • • S Ns ) ks ) (37) 

and the four functions = (3 kr , ^(1) = ia kr , ^(-1) = a kr , and = -ip kr , 

compare eq. (0). The eigenstates for fixed number of charge excitations Nc can then be 
written as 

\ki, n; . . . k Nc ,r Nc ; (Ci, . . . Cjv )fc ; (Si, . . . SN s )k s ) = 
~n c 



r=l 



|fci, . . . fcjv c ; (Ci, . . . Cjsr c )k a ] (Si, . . . SN 3 )k s ) ■ (38) 



This state corresponds to (Nc + J2 r T r) [{Nc — Z)r r r)/2) spinless Fermions in the upper 
(lower) Peierls subband. 

16 



The corresponding energies and momenta of these states are obviously given by 

N c 

E = Y1 E(kj + $cs)Tj + UN d (39a) 
3=1 

N C 

P = k s + Y,(kj + ®cs) modvr/a (39b) 

3=1 

where E(k) has been given in eq. (H). 



C. Band picture interpretation of the spectrum 

1. Translational invariant case 

The exact solution for the Harris-Lange model can be interpreted in terms of upper and 
lower Hubbard bands. To simplify the discussion we will ignore the momentum shift §cs m 
this subsection. For U > W the ground state of the half-filled band N = L has energy zero 
and is 2 i -fold degenerate, and we may choose the fully polarized ferromagnetic state as our 
reference state, |FM) = | f, . . . f). Note that this state has momentum ir/a on a chain with 
an even number of sites L, see eq. (|35b| ). 

We may now add an electron. We obtain all exact eigenstates for N = L + 1 electrons, 
Nd = 1, and all spins up as 

\k; (.) fec=0 , (T, . . . T) fcs=0 ) = 4J FM > (40a) 



with momentum P = k + ir/a and energy E = e(k) + U. The state in eq. ( 40a|) is inter- 



preted as a particle at momentum k in the upper Hubbard band which has the dispersion 
relation e(k) + U. The momentum ir/a is attributed to the ferromagnetic reference state. 

We may equally well take out an electron from the fully polarized state. We obtain all 
exact eigenstates for N = L — 1 electrons, = 1, and all spins up as 

\k; (o) fcc=0 , (t, . . . T) fcs =o) = -c^ /a _ fclT |FM) (40b) 

with momentum P = —(it /a — k) + it /a and energy E = e(k). Note that the states of 
eq. ( f40b|) and those of eq. ( |40a|) can be generated from each other by the particle-hole 
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transformation of eq. ([161) . Their momenta differ by it /a since their respective numbers of 
double occupancies differ by one. 

Since the ground state corresponds to a completely filled lower Hubbard band we interpret 



momentum it /a is again attributed to the ferromagnetic reference state. The lower Hubbard 
band must have the dispersion relation e(k) for particles because a hole at kh = tt/cl — k has 
momentum P = —(n/a — k) and energy E = —e(kh) = —e(—k + it /a) = e(k). 

The band structure for the Harris-Lange model is depicted in figure 0. It displays the 
parallel upper and lower Hubbard bands with band width W separated by a distance U. 
It is amusing that the celebrated Hubbard-I approximation (Hubbard 1963); (Mazumdar 
and Soos 1981) also gives parallel bands. Those bands, however, carry a spin index while 
charge-spin separation is most essential in one dimension. Furthermore, the width of those 
bands is only half of the exact band width W. 

Our band structure picture has to be used carefully if there are more than one double 
occupancy or hole. Figure [I] suggests that there are L states available both in the upper and 
in the lower Hubbard band, altogether 2L independent states. However, this cannot be the 



appropriate exclusion principle between particles in the upper Hubbard band and holes in 
the lower Hubbard band can be formulated. 

For fixed spin background and fixed hi, k 2 there are four exact eigenstates with two 
charges at k\ ^ k 2 . They all have the kinetic energy T = e(ki) + e(k 2 ). They correspond 
to four different charge excitations in the band picture: (i) two particles at momenta k%, 
k 2 in the upper Hubbard band, (ii) two holes at momenta irja — k±, it /a — k 2 in the lower 
Hubbard band, (iii) a particle at momentum k\ in the upper Hubbard band and a hole at 
momentum irja — k 2 in the lower Hubbard band, (iv) a particle at momentum k 2 in the 
upper Hubbard band and a hole at momentum irja — kx in the lower Hubbard band. The 
condition k\ ^ k 2 is naturally fulfilled in cases (i) and (ii), if we assign a fermionic character 



the state in eq. ( |40b| ) as a hole in the lower Hubbard band at kh 



it /a — k, and the 
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to our particles in the upper and holes in the lower Hubbard band, respectively. In case (iii), 
however, we have to explicitly demand that the momentum at which we create the hole, 
kh = Tr/a — k%, fulfills k\ 7^ &2, i- e., kh 7^ 7r / 'a — k\. This is the same condition which results 
from case (iv). 

We thus see that a particle in the upper Hubbard band at momentum k actually blocks 
the creation of a hole in the lower Hubbard band at momentum n/a — k (this is probably the 
simplest example of a "statistical interaction" , see (Haldane 1991)). With this additional rule 
the counting of states in the band picture is correct, and the band picture interpretation 
gives indeed the exact results for the Harris-Lange model. The effective Hamiltonian for 
fermionic particles in the upper (u k ) and lower (l k ) Hubbard band thus reads 

H^ d = P ul E [(U + e(k))nt + e(k)h l k ] P ul (41a) 

\k\<n/a 

P ul= II I 1 - (l - K/a-k) "*] ( 41b ) 
\k\<Tr/a 

with hi = u k ii k , n l k = l k l k . The projection operators guarantee that there is no hole in the 
lower Hubbard band at momentum n/a — k, if there is already a particle at momentum k 
in the upper Hubbard band. For half-filling at zero temperature the lower Hubbard band is 
completely filled. 

2. Dimerized Harris-Lange model 

The case of the dimerized Hubbard model can be treated accordingly. The upper and 
lower Hubbard band now split into two Peierls subbands with dispersion relations ±E(k). 
Formally the band structure Hamiltonian for the lower band becomes (compare eq. (0)) 

f^(6) = £ e(k) (f+f fc - it +7T/ J k+7T/a ) - iA(k) (j+ +7r/ J k - i+i k+n/a ) , (42) 

|Jfc|<7r/(2o) 

and a similar expression holds for the upper Hubbard band. 

The band picture Hamiltonian can easily be brought into diagonal form as in the Peierls 
case. The quasi-particles in the four subbands are finally described by 
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< m ' band = P u+ i + Pn-i- E (U + E(k))nl + + (U-E(k))nl_ 

|A|<7r/(2o) 



+ E(k)n{, - E{k)n l u_ P u+l+ P u -i- (43) 



iW = n 

|fe|<7r/(2a) 



1 " "U, ± ) nJU 



with ± = ut±u k ± , n' fe ± = Z^" ±l k ± as the number operators for the quasi-particles for the 
upper (r = +) and lower (r = — ) Peierls subband in the upper (u) and lower (Z) Hubbard 
band. The quasi-particles in each subband obey a fermionic exclusion principle in the same 
Hubbard band. In addition, a particle at momentum k in the upper Hubbard band in the 
upper (lower) Peierls subband blocks the creation of a hole at momentum — k in the lower 
Hubbard band in the upper (lower) Peierls subband. There is no hole in the lower Hubbard 
band at momentum — k in the upper or lower Peierls subband, if there is already a particle 
at momentum k in the upper Hubbard band in the corresponding Peierls subband. Note 
that the reciprocal lattice vector is now given by it /a. Thereby the projection operators 
guarantee the proper counting of states. 

The resulting band structure is shown in figure The upper and lower Hubbard band 
are both Peierls split and display the Peierls gap W5 at the zone boundaries ±7r/ (2a). Note 
that the upper (lower) Peierls subbands are still parallel. 

D. Band picture interpretation of the current operator 

1. Translational invariant case 
According to the spectral decomposition of the current- current correlation function, 



eq. (|22|), we need to determine the excitation energy E n — E of an exact eigenstate \n) 
and its oscillator strength |(0|j|n)| 2 . The respective total momenta of these states are P 
and P n . 

We are interested in optical excitations from a state with singly occupied sites only. 
The excited states which can be reached from this state have one pair of hole and double 
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occupancy, 1. e., 

\0) = \(S 1 ,...S L ) ks ) (44a) 

|n) = \h, k 2 ; (mo) kc=0 ; . . . S L - 2 )k0 (44b) 

where we used the fact that j creates a charge triplet with kc = 0. Note that (k\,k 2 ) is the 
same state as (k 2 , k\). We denote k\ = k + q/2, k 2 = ft/a — k + g/2 since we will finally 
represent the state \n) by a particle in the upper Hubbard band at momentum k + q/2 and 
a hole in the lower Hubbard band at momentum k — q/2. 

Since j conserves the total momentum we already know that Pq = P n which implies 
ks = q + k' s (L — 2)/L, see eq. ( |35b|) . Hence, k' s = L{ks — q)/(L — 2) has to hold. Recall 
that k' s is quantized in units of 27r/((L — 2)a). These considerations imply that the charge 
(spin) sector in |n) carries momentum q (—q) relative to |0). In the thermodynamical limit 
the excitation energy is given by 

E(k, q) = U + e(k + q/2) - e(k - q/2) = U + At sin(ka) sin(ga/2) . (45) 

Note that the excitation energy does not depend on the spin configuration. For this reason 
it is possible to find a formally equivalent band picture for the charge sector alone that gives 
the same optical absorption as the original model. Since j itself carries all the information 
on the conservation laws (momentum, charge, and spin quantum numbers) we may equally 
well work with the (normalized) states 

|£ _|_ S.-- - k + -) = — f e i(k+q/2)l 1 a e i(n/a-k+q/2)l 2 a _ e i(k+q/2)l 2 a ^(tt / a~k+q/2)l 1 a\ 

(-l) h \S[, . . . S'^, m hi S' h , . . . S' l2 _ 2 , o l2 , <4_ 1; . . . S' L ^ 2 ) (46) 



rather than the exact eigenstates of eq. (p4f) . This will simplify the notation since we do not 
have to take any summation restrictions into account. 

For fixed (k, q) and fixed spin configuration (S[, . . . S' L _ 2 ) we calculate 

(0&..-I* + %\ ~ k + |) = -teae(k)e^ 2 (47) 

j E e^ la (0\S[, . . . SU (Hi+i - Uh+i) , s' h . . . S' L _ 2 ) . 
L I 
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We define the operators x+ and x via their product 



xZx 



E 71 E e^- r)a (0|^, . . . 5j_ lf (Tdw - Uh + i) , • • • S' L _ 2 ) 



5(,...5^_ 2 ^ i,r 



(48) 



(51 



L-11 



S' r i (Ir+ltr — tr+lir) > 5 r _ 1 , . . . S[\0) 



Summed over all intermediate spin configurations the oscillator strength for fixed (k, q) now 
becomes 



9, 



—ieae(k) 



X q X q 



(49) 



It is clear that we have hidden a very difficult many-body problem in the operators x g . 

Nevertheless we are now in the position to identify the interband current operator in the 
band picture. It is given by 



~band 
Jinter 



E ^eae(&0 {^k+q^k-q^q K-q/2^k+q/2 X q 
\k\,\q\<n/a 



(50) 



This operator acts in the same space as the band Hamiltonian of section [111 (J| . It is seen 
that the condition k + q/2 7^ ix/a — (k — q/2) is automatically fulfilled since e(7r/(2a)) = 0. 
Consequently, the projection operators in eq. (|4lD can again be ignored for the case of linear 
optical absorption. 



2. Dimerized Harris-Lange model 

For a distorted lattice the current operator can also modify the momentum of the state 
by ix /a. Thus we address four possible states for fixed (k,q) from the reduced Brillouin 
zone, \k + q/2; ix/a - k + q/2), \k + q/2; —k + q/2), \ix/a + k + q/2; ix/a - k + q/2), and 
\ix/a + k + q/2; — k + q/2). The same analysis as in the previous subsection leads us to the 
definition of the operators x^(6, rj), x q (5, rf) with the property 
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x+{6, v)x qf (S, V )= £ -1 £ e^-^ a (l + ^ + + 77))(1 + ^ + + r?)) 

s' x ,...s' L _ 2 lj l,r 

(o\s[,... su (Ui+i - Uh+i) , si, . . . s' L _ 2 ) (si) 

(S'l-21 ■ ■ ■ S' r , (jr+ltr ~ Tr+llr) > 'S'r-l' • • • 5*1 1 ) . 

In practice, q' = q or q' = q + n/a. 

The interband current operator becomes 

"band "band _i_ "band 

Jintcr Jintcr,+ ' Winter,— 

"band ( "band \ ~'~ 

Jinter,+ \Jmter,— J 

Jinter,- = E | _ ^ae(k) lt- q /2^k+q/2 ~ ^k-q/2+TT /a^k+q/2+n /a %q (52) 

[g|,|fe|<7r/(2a) k 



A(Jfe) 

+ ea 



^k-q/2^k+q/2+Tr/a ^k-q/2+ir / 'a^k+q/2 



X q+Tr/a 



5 

Again the conditions k + q/2 ^ -n j a — ik — qjl) and k + q/2 ^ — {k — q/2) will not be violated 
since e(7r/(2a)) = and A(0) = 0, respectively. Consequently, the projection operators in 
eq. ([43]) can be ignored for the case of linear optical absorption. 

In semiconductor physics one prefers to work with the dipole operator rather than the 
current operator to set up the perturbation theory in the electrical field (Haug et al. 1990). 
The corresponding expressions for the dipole operator for Hubbard interband transitions are 
derived in appendix V 



IV. OPTICAL ABSORPTION IN THE HARRIS-LANGE MODEL 

A. Spin average 

All states with no double occupancy are possible ground states in the Harris-Lange model 
at half-filling. Instead of looking at the optical absorption for a specific state |0) it is more 
reasonable to calculate the average absorption, i. e., 
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Im{xH} = ^53 Im {X|o>H} • 



(53) 



|o> 



For the Hubbard model this corresponds to temperatures k B T ^> J = 0{W 2 /U) ("hot-spin 
case"). The calculation is performed in appendix |A5|. We find the result 



1 

2L 



Sq,q> 



:i + Sr l ) 2 g(q) + (S + r l ) 2 g(q + 



7T s 



+ S m > + K(l + 8rj)(S + rj) 
3 



7T , 



g(q) + g(q + -) 

a 



(54a) 
(54b) 



5 + 4 cos(ga) 

The spin problem could thus be traced out completely. It is seen that x q keeps its operator 
character until we have expressed the current operator in terms of the Fermion operators 
for the four Peierls subbands. 



B. Optical absorption 

1. Translational invariant case 
The real part of the average optical conductivity becomes 

Re{a(co > 0)} = £f±- £ (eae(k)) 2 g(q)5(uj - E(k, q)) (55) 

|g|,|fc|<7r/o 

with E(k,q) from equation (|4"5|). The above expression can be written as 

t r 1 f 1 dx \lx 2 — u 2 3 

^ >0) = ^ i, ^Vl-x 2 9-8x' (56) 

for the reduced optical conductivity with u = \uj — U\/W < 1. This integral can be expressed 
as a sum over elliptic integrals but we rather prefer to discuss some special cases. 

The optical absorption is restricted to \uj — U\ < W. Near the band edges the absorption 
increases linearly which can be seen from equation (|5^) by a transformation x 1 — y/\u\. 
The more interesting case is ou — U. Now the integrand displays a 1/x singularity for \u\ — > 0. 
The parallel Hubbard bands give rise to a logarithmic divergence, a{u — > U) ~ | \n(u — U)\, 
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since their large joint density of states for u = U survives even in the presence of a spinon 
bath that provides any momentum to the charge sector. 

The overall behavior of the optical absorption is shown in figure 0. The same absorption 
curve has been obtained earlier in (Lyo et al. 1977) for their "random" spin background. 
The result for their "ferromagnetic" spin background follows when we put g(q) = 1, as 
expected. We will comment on their Neel state results in (Gebhard et al. Ill 1996). 



2. Dimerized Harris-Lange model 



We have to diagonalize the interband current operator of eq. fl52|) in terms of the Peierls 
operators for the lower Hubbard band 



h — o>kh,- + fikh 



h+w/a — —ifikh- + iockh,^ 



(57a) 
(57b) 



for \k\ < 7r/ (2a). The transformation for the upper Hubbard band is analogous. With this 
definition the Hamiltonian in the band picture interpretation became diagonal, see eq. ([43D. 
The interband current operator becomes 



"band 
Jinter,- 



E E X T,r>{k,q)i+_ q/2jT u k+q/2y 

T,T>=±l\k\,\q\<ir/(2a) 



(58) 



with 



\ +i+ (k, q) = iea 
A +> _(&;, q) = iea 



-e(k)(a + (3*_ + p + al)$* + ^(a + a*_ - (3 + (3*_)x^ /a 



(59a) 
(59b) 



and X--(k, q) = — X +t+ (k, q), X- 7+ (k, q) = X + -(k, q). Here we used the short-hand notation 
a± = cik± q /2 etc. Note that these quantities can be complex for q ^ 0. 
The average optical conductivity becomes 



Re{cx(co> > 0,5,77)} 



Lau 



E \X T y(k,q)\ 2 6(uj-E T y(k,q)) (60a) 



t,t'=±1 | 9 |,|fe|<7r/(2a) 
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with the absorption energies between the respective Peierls subbands 



E T y (k, q) = U + r'E{k + q/2) - rE{k - q/2) , 
see eq. (f|3]) and figure ||. The transition matrix elements are given by 



(60b) 



\K,r'{k,q)[ 



2L 



9{q) 



'l + 5 v )e(k)f T)T/ + TT'(5 + r])^p-f Tr 



+ g(q + 



7T . 



A(k) 

(5 + V )e(k)f T y + tt'{1 + Sri)-y-f T> ^ 



(60c) 



with the help functions 



/+,+ (*> ?) = /-,-(*>?) = a k+ g /2 a l- g /2 ~ Pk+q/2Pk-q/2 
f+ -0> Q) = /-,+ (*> g) = ^k+ q /2Pl-q/2 + Pk+ q /2^t- q /2 



(61a) 
(61b) 



where (3k are given in eq. (^). It can easily be checked that the case 5 = r\ = is 
reproduced. For 5 = 1, r] — one recovers the result for the average optical conductivity 
of L/2 independent two-site systems since E(k) = 2t, X + -(k,q) = 0, and \X +t+ (k,q)\ 2 = 
(2tea) 2 (g(q)+g(q + 7c/a))/(2L): 



Re{a(u>0,5 = l,r ] = 0)} = ^7r5(u - U)^^ 



(62) 



2 Lau " y " ~ ' 4 

where we used J*^ dqj (2-k) g(qa) = 1. For the direct calculation we have to recall that only 
the singlet of the four spin states contributes, and the hopping between the two sites is 2t. 

For general S, r] it is necessary to evaluate the optical conductivity in eq. (^) numerically. 
An example is shown in figure [|. It is seen that now there are two side-peaks in the optical 
absorption spectrum. The new peaks due to the Peierls distortion are weaker than the one 
at oj = U and vanish for S — > 1. For large lattice distortions the dominant contribution to 
the side peaks comes from the small-g transitions between different Peierls subbands. Their 
oscillator strength is maximum for u = U±W\/S which determines the position of the peaks 
for large 5. The Peierls gap between the bands shows up in the optical spectrum. For small 
lattice distortions all (k, q) contribute. The signature of the Peierls gap is smeared out and 
the position of the side peaks cannot be expressed in terms of a simple function of 5. 
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V. OPTICAL ABSORPTION IN THE EXTENDED HARRIS-LANGE MODEL 



A. Extended dimerized Harris-Lange model 

Strongly isotropic, almost ideal one- dimensional systems like polymers or charge-transfer 
salts are not properly described by the Hubbard model of eq. for two reasons: (i) the 
Peierls distortion is not taken into account, and (ii) the residual Coulomb interaction between 
charges beyond the Hubbard on-site interaction is neglected. The exponential decay of the 
Wannier wave functions naturally allows to limit the interactions to on-site and nearest- 
neighbor Hubbard terms. In the "Zero Differential Overlap Approximation" it is further 
assumed that only the direct Coulomb term has to be taken into account for the nearest- 
neighbor Coulomb interaction (Kivelson, Su, Schrieffer, and Heeger 1987); (Wu, Sun, and 
Nasu 1987); (Baeriswyl, Horsch, and Maki 1988); (Gammel and Campbell 1988); (Kivelson, 
Su, Schrieffer, and Heeger 1988); (Campbell et al. 1988); (Painelli and Girlando 1988); 
(Painelli and Girlando 1989); (Campbell, Gammel, and Loh 1990). 

Optical absorption spectra for the extended dimerized Hubbard model could only be cal- 
culated numerically for small system sizes. Within such an approach the Hamiltonian is ex- 
plicitly diagonalized (Soos and Ramesesha 1984); (Tavan and Schulten 1986); (Guo, Mazum- 
dar, Dixit, Kajzar, Jarka, Kawabe, and Peyghambarian 1993); (Guo, Guo, and Mazumdar 
1994). Since the dimension of the Hilbert space increases exponentially (dimff = 4 L ) the 
numerical analysis is restricted to short chains (L < 12) due to the limited computer power. 

Therefore, it is natural to analytically investigate the (dimerized) Harris-Lange model 
with an additional nearest-neighbor interaction. We will show below that the optical spec- 
trum can still be calculated for this model which is equivalent to the extended dimerized 
Hubbard model to order t(t/U)~ x , t(t/U)°, and t(V/U)°. The dimerized extended Harris- 
Lange model reads 



27 



V = Y,(ni-l)(ni +1 -l). (63b) 
i 

For half-filling the ground state of the extended dimerized Harris-Lange model is still 2 L -fold 
spin degenerate because every site is singly occupied for \V\ < U/2. The energy of these 
states is zero, E = 0, irrespective of the dimerization value 5. 

The double occupancy and the hole in the excited states now experience a nearest- 
neighbor attraction while the spin sector remains unchanged. Thus we may immediately 
translate V into our band picture as 

^ band = -f E ™s(qa) £ u + k+q u k lJU (64) 

\q\<n/a \k\,\p\< n / a 

which describes the scattering of a hole in the lower Hubbard band with a particle in the 
upper Hubbard band. Again, the projection operators can be disregarded for the optical 
absorption. 

B. Equation of motion technique 

In the presence of the nearest-neighbor interaction it becomes increasingly tedious to 
separately calculate the exact eigenenergies and oscillator strengths. We rather prefer to 
directly calculate the optical conductivity from an equation of motion approach. 

1. Translational invariant case 

Since we are interested in the real part of the optical conductivity we can concentrate 
on the particle current density, 

a> = ^«0(*)lil*(*)>+h.c.) (65) 

where |0(t)) and \^{t)) are the time evolution of the ground state with and without the 
external perturbation and thus obey the corresponding Schrodinger equations. We have 
already used the fact that we want to calculate the linear absorption. 
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We write 



(o(Ob«0> = E-^(^)^(°Wlt,/2% + ,/2l*W) 

= E A (MW*)- (66) 

Upon Fourier transformation we obtain 

&> = ^ (e A ( fc > + • (67) 

\k,q ) 

Since we are interested in the optical conductivity for positive frequencies (optical absorp- 
tion) we may disregard the second term which contributes to u < 0. 
The equation of motion for jk;q(t) becomes 

l dj^t)_ = {m [l+_ q/2 u k+q/2 ,H* T «]_ m) ) - ^(0(t)|t g/ A +9/2 j|0(t)) (68a) 
uj kiq (u) = E(k, q)j k;q (u) - 2V (cos(ka)j c q (u) + sm(ka)f q (cu)) - A+(fc, q)^ (68b) 

where we kept the expansion linear in the external perturbation, and performed the Fourier 
transformation. Furthermore, we introduced the abbreviations 

w-isKS)*^- (69) 

For our calculations we only need j q {oS) since our current operator preserves parity. The 
particle current density for positive frequencies becomes 

CL>o) = — (2tiea) E * + q j» ( 7 °) 
a q 

which is proportional to the external field. 

We introduce the function 

_ 2 (cos kaf 

F ^-L^^m¥) (71) 

which allows us to finally express the optical conductivity as 

The result will be discussed in the next subsection. 
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2. Extended dimerized Harris-Lange model 



The same procedure can be applied to the dimerized case where it is best to start from 
the diagonalized Hamiltonian in the form of eq. (fH|), and the current operator in the form 
of eq. (|58p . The calculations are outlined in appendix |B[ 

We introduce the three functions F 123 as 



cos 2 (/ca) 



|fc|<7r/(2o) 



F 2 (q) 



sin 2 (fca) 



|fc|<7r/(2o) 



1 



+ 



l 



u — E-- uj — 

1 1 

+ 



u - E- 
1 



+ 



uj — E +i 
1 \ 



LJ - E +,+ , 



+ \fi 



1 



+ 



F 3 (q) 



L 



cos(ka) sm(ka)l /++/+ 



uj — £L.+ uj — 
1 



|fc|<7r/(2a) 



uj — E j 



/;.+/+.- 



uj- E_ 



+ 



uj-E,_ 



(73a) 



(73b) 



(73c) 



where / TjT / = f T y(k,q) and i£ T)T / = E TT i(k,q) were introduced in eq. (pT|) and eq. 
Furthermore, we abbreviate A; = (1 + 8r))Fj — (77 + 5)F 3 , = (5 + r^)Fj — (1 + 5r])F 3 
(j = 1, 2), d = (1 + ^) 2 F! + (5 + r]) 2 F 2 - 2(1 + <fy)(<S + 77)^3, and C 2 = (1 + 5r/) 2 F 2 + (5 + 
F\ — 2(1 + 5r/)(<5 + r/)^. The real part of the average optical conductivity can then be 
expressed as 



Re{a(uj > 0, V, 5,r])} = Re{a(uj > 0, 5, 77)} 

WV ± (WW 2 



16au;L 



to £71 



1 



|<z]<7r/(2a) 



VFOCl + yFa)-^) 1 



^(g) A 2 + E 2 + y (F 1 F 2 - Fi)C x 



71 , 



(74) 



+ g(q + -) [A 2 + Bf + V(F 1 F 2 - F 2 )C 2 
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The result for V = is given in eq. fl60|). In the following we will discuss the results for the 
average optical absorption in the presence of a nearest-neighbor interaction. 



C. Optical absorption 

1. Translational invariant case 

The help function F(q) can be calculated analytically with the help of eqs. (2.267,1), 
(2.266), and (2.261) of (Gradshteyhn and Ryzhik 1980). The result is 

u-U - ^f(uj - U) 2 - (4tsm(qa/2)f for \u - U\ > |4tsin(ga/2)| 
[4tsm(ga/2)] 2 | sin ( ga / 2 )) 2 _ ( w _ u)z f or \ u - U\ < |4tsin(ga/2)| 

(75) 



F(q) 2 



The result for V = 0, eq. (|5q) , follows after the substitution of x — sin(ga/2) into equa- 
tion (0). 

The total optical absorption is shown in figure |5|. For arbitrarily small V > there is 
a bound exciton which is the standard situation for one-dimensional short-range attractive 
potentials between a positive (hole) and negative charge (double occupancy). This is evident 
from the form of F(q) which allows for excitons with momenta q if \u — U\ > \4t sin(ga/2)| 
which is fulfilled for q = for all V > 0. 

When the attraction between the two opposite charges is strong, the full exciton band 
with width W exc = 4t 2 /V is formed. This can be seen from the zeros of the denominator in 
eq. © in the region \u - (U — V - W cxc /2)\ < W cxc /2. One finds from 1 + VF(q) = that 

to = U — V 1 — cos qa . (76) 

This is precisely the dispersion relation for bound pairs in one dimension with nearest- 
neighbor hopping of strength t exc = t(t/V): at large V the excitons are essentially nearest- 
neighbor pairs of opposite charges which coherently move with the hopping amplitude t exc . 
Note that this motion requires an intermediate ("virtual") configuration where the two 
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charges are not nearest neighbors. Consequently, the hopping integral of the individual 
constituents, t, is reduced by the factor t/V for the motion of the pair. The full band 
becomes apparent when W exc + V > W or V > W/2, see figure |5]. Recall that the spin 
sector provides any momentum to the charge sector. The momentum transfer, however, is 
modulated by the function g(q) which is maximum at q = it /a and reflects the fact that 
states with antiferromagnetic spin correlations are best suited for optical absorptions since 
they contain many neighboring singlet pairs. Hence, the q = 7r/a-exciton dominates over 
the q = 0-exciton for V > W/2. 

It is amusing to see that the optical absorption of a Peierls insulator and a Mott-insulator 
(extended Harris-Lange model: U ^> W, V > W/2, J = 0) can look very similar, compare 
figure 2 of I and figure [|. This has already been noted long time ago by Simpson (Simpson 
1951); (Simpson 1955); (Salem 1966); (Fave 1992) who explained the optical absorption 
spectra of short polyenes in the above exciton model. It is seen that Simpson's model is 
naturally included in our strong- correlation approach. For real polymers, however, Simpson's 
original approach is not satisfactory. A fully developed exciton band only exists in the 
presence of an incoherent spin background. Now that even the spin- Peierls effect is excluded 
one can by no means explain the Peierls distortion of the lattice as an electronic effect. This 
does not exclude other, e.g., extrinsic, explanations for a lattice distortion. 



2. Extended dimerized Harris-Lange model 

The full spectrum has to be determined numerically. An example for various values of V/ t 
is shown in figure |6]. For large V/t we have a fully developed exciton band which is itself 
Peierls-split into two branches. Thus one obtains four van-Hove singularities in the optical 
absorption spectrum. The Peierls gap is given by A^ xc = 6W exc = 4t 2 5/V. Even for V = W 
it is smeared out since V/t is not too large yet and the phenomenological damping 7 is 
already of the order of the gap. 

For small V/t we obtain the signature of the q = exciton for 5 = 0, V > 0. For 
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intermediate V this peak develops into a van- Hove singularity of the upper exciton subband. 
The signatures of the second van-Hove singularity of the upper band are clearly visible 
for V = W/2. The peaks of the Peierls subbands for V = 0, 5 ^ are both red-shifted. 
The peak at lower energy increases in intensity and finally forms the lower exciton subband 
while the peak at higher energy quickly looses its oscillator strength. 

VI. SUMMARY AND OUTLOOK 

In this paper we addressed the optical absorption of the half-filled Harris-Lange model 
which is equivalent to the Hubbard model at strong correlations and temperatures large 
compared to the spin energy scale. It is extremely difficult to analytically calculate optical 
properties of interacting electrons in one dimension. For strong coupling when the Hubbard 
interaction is large compared to the band-width matters considerably simplified since the 
energy scales for the charge and spin excitations are well separated. We were able to derive 
an exactly equivalent band structure picture for the charge degrees of freedom and found 
that the upper and lower Hubbard band are actually parallel bands with the band structure 
of free Fermions. We have taken special care of the spin background which can act as a 
momentum reservoir for the charge system. Since we can exactly integrate out the spin 
degrees of freedom for the Harris-Lange model we were able to solve the problem even in the 
presence of a lattice dimerization and a nearest-neighbor interaction between the electrons. 

For a vanishing nearest-neighbor interaction we found a prominent absorption peak 
at uj — U and additional side peaks in the absorption bands \uj — U\ < W in the pres- 
ence of a lattice distortion 5. When a further nearest-neighbor interaction between the 
charges was included, we found the formation of Simpson's exciton band of band-width 
W cxc . = At 2 /V for V > W/2 which is eventually Peierls-split. As usual the excitons draw 
almost all oscillator strength from the band. 

It should be clear that the Harris-Lange model with its highly degenerate ground state is 
not a suitable model for the study of real materials. The results presented here are relevant 
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to systems for which the temperature is much larger than the spin exchange energy. Real 
experiments are not carried out in this "hot-spin" regime but at much lower temperatures 
for which the system is in an unique ground state with antiferromagnetic correlations. Un- 
fortunately, this problem cannot be solved analytically. In the third and last paper of this 
series (Gebhard et al. Ill 1996) we will employ the analogy to an ordinary semiconductor 
(electrons and holes in a phonon bath) to design a "no-recoil" approximation for the char- 
geons in a spinon bath. It will allows us to determine the coherent absorption features of 
the Hubbard model at large U/t. 
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APPENDIX A: THE HARRIS-LANGE MODEL 



1. Sum rules 



We briefly account for the sum rules. We have 

Mi 



roo A/ . o Kf 

/ dulm{ X (uj)} = n-^Y: <°b» =TT i W|0>. 
Jo La „ 1 La 

It is a standard exercise to show that 

(0|j 2 |0) = (2tea) 2 £ (l + (-1)^)" (l + (-l)^) 2 (0| (~ - S,S m ) |0) . 



We define the positive quantities 



C 



even, odd 



^^(oifi-s,s i+l) io> 



4 



and may then write 

roo 

/ duj Im{xM} = nJ\f ± a(2tef 
Jo 



1 + 5f (1 + r]f C| ven + (1 - <5)" (1 - t]) 2 C° s dd 



(Al) 



(A2) 



(A3) 



(A4) 



If we average over all possible states |0) we obtain (7° ven ' odd = \j% since only the singlet 



configuration contributes. Hence, 



du\m{x{uj)} = irjV ± a(teY (1 + 6rjY + (5 + rj) 



The area under the curves for a re d(uj), eq. fl2"4T) , are thus given by 



(A5) 



o ^ rcdM = 4 



1 + 5)" (1 + r]f C° ven + (1 - 5f (1 - ??) 2 Cg dd 



uF ff «dM = ^ (l + Srjf + iS + rjy 
o W lb 1 



(A6a) 
(A6b) 



2. Momentum of eigenstates 

First we will assume that the number of sites and the number of particles is even. 
Nc and N$ will thus also be even. Let T be the translation operator by one site, i. e., 
Tc/ )Cr T _1 = ci + x a . We have to show that 
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f |*> = e"* Pa |^) (A7) 

holds. We have to distinguish two cases: (i) a spin is at site L in |\&) or (ii) a charge is at 
site L in |\&). 

a. case (i): the operator T shifts the spin from site L to the first site. This results in a 
phase factor (—1) since one has to commute the Fermion operator (N — l)-times to obtain 
the proper order in Furthermore, the states in |\l>) have, relative to those in T|\&), 

1. shifted the spin sequence by one unit. This results in a phase factor exp(i/csa); 

2. an additional factor (—1) for each doubly occupied site. This gives a phase fac- 
tor exp(i7!-Nd); 

3. a Slater determinant in which each site index is shifted by one. This results in a phase 
factor exp(i T,f=i(kj + $cs)a)- 

In sum we obtain 

_ e -iPa e ik s a e iTrN de iY^= 1 ( k 3+ iS 'Cs)a 

N c 

P = k s + (7r/a)(N d -l) + Y,(k j + <S> cs ) mod 2ir/a . (A8) 

j'=i 

6. case (m,): the operator T shifts the charge from site L to the first site. The states 
in |\&) have, relative to those in T|^), 

1. shifted the charge sequence by one unit. This results in a phase factor exp(ik c a); 

2. an additional factor (—1) for each doubly occupied site. This gives a phase fac- 
tor exp(irrNd); 

3. a Slater determinant in which 

(a) each site index is shifted by one. This results in a phase factor exp(iJ2f=i{kj + 
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(b) the last row and the first row are interchanged. This gives an additional factor 
(_l)^c-l = _ 1; 

(c) the first row is (1, . . . 1) instead of (exp(i(ki + §cs)La), . . . exp(i(k Nc + § C s)La)) 
= (1, ... 1) exp(i(k c — k c )a). This gives an additional phase factor exp(— i(k c ~ 
k s )a). 

In sum we obtain 

I _ e -iPa e ikca e inN de iJ2^=i( k J+^Cs)a.^_^ e -i(kc-k s )a 

N C 

P = k s +(Tr/a)(N d -l)+Y,(k j + ®cs) mod Ivja (A9) 

as before. If the number of particles N is odd, the momentum is shifted by another factor 
of n/a, if one repeats the above arguments. This proves that the total momentum of the 



state is indeed given by eq. (|35b|). 



3. Energy of eigenstates 

We want to prove that 

N c 

H nh \^) = E\m) ; E = ^2e(k j + $cs) + UN d . (A10) 

Again we restrict ourselves to even N. The bulk terms are simple since there is no hopping 
across the boundary. A hopping process of a double occupancy and a hole are equivalent: 

(-l)'l . . . a l+1 ...)>-> -(-1)1 • • • °i •/+! • • •) = (~1) ,+1 | ...d • • •) (Alia) 
| . . . oj a i+ i . . .} i->- | . . . ai o l+1 ...) . (Allb) 

The extra minus sign which appears when a double occupancy moves has been taken care 
of in the wave function by the phase factor (—1)' for a double occupancy at site I. Now 
that there is no difference in the motion of double occupied sites and holes they dynamically 
behave as spinless Fermions. Since the Slater determinant is the proper phase factor for 
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non-interacting Fermions eq. ([A10Q holds for the bulk terms. This also shows that only 
the Harris-Lange model with hopping amplitudes |£lhb| — I^uhbI can be solved. Another 
integrable but trivial case is tLHB = 0, £uhb 7^ and vice versa. 

We now address the boundary terms. A typical configuration for which transport across 
the boundary is possible is 1 5*1, . . . C£). The phase of the configuration is given by a Slater de- 
terminant in which the last row has the entry (exp(i(fci + $cs)La), . . . exp(i(kN c + §cs)La)) 
= (1, ... 1) exp(i(kc — ks)La). The action of Hhl moves the charge from site L to the first 
position by which an extra minus sign occurs since the electron operator for the spin had to 
be commuted with (N — 1) other electron operators. These phase factors have to be com- 
pared to the corresponding configuration in E\ty). Relative to the configuration in Hul\^) 
it has 

1. shifted the spin sequence by minus one unit. This results in a phase factor exp(— ikga); 

2. shifted the charge sequence by one unit. This results in a phase factor exp(ikca)] 

3. a Slater determinant which has (exp(i(fci + <&cs) a )-> ■ ■ ■ exp(z(fc A r c + $cs)°0) i n the first 
row. 

The boundary terms should give the same result as the bulk terms. This leads to the 
condition 

— e i( . k C- k s)La _ e -iksa e ik c a{^N c -l (A12) 

which is obviously fulfilled. The proof for odd iV is analogous, and eq. (|A10 ) holds for all N. 



4. Electrical dipole operator for the Harris-Lange model 

We can derive the electrical dipole operator for the Harris-Lange model from its definition 
in equation (A. 22) of I. We use the Hamilton operator in the band picture interpretation, 
eq. (|4"T|), and the corresponding current operator, eq. fl50|). Since x q can be replaced by its 



average value ^ g(q) / (2L) one easily sees that the dipole operator becomes 
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inlcr E /^interC^' Q) ( 

\k\,\q\<n/a 



U k+q/2^k-q/2 + ^k-q/2^k+q/2 



Winter (^5 <?) ~ * 



A(fc,g) 



(Al3a) 



(A13b) 



#0,g) V 2L 

with E(k, q) = U + e(/c + q/2) — e(k — q/2). One sees that the dipole matrix element is of 
the order t/U as it should be for interband transitions. Furthermore, for small momentum 
transfer we obtain 



ACUfc;?->o)~e(fc). 



(A14) 



This is the correct form since we create a neighboring hole and double occupancy which 
corresponds to an electric dipole between nearest neighbors. 

The procedure is readily generalized to the dimerized Harris-Lange model. The interband 
current operator in terms of the Fermion operators for the four Peierls subbands is given in 
eq. (pBj), and the diagonalized Hamiltonian can be found in eq. ((43|). One readily finds 



jpdim. HL 
inter 



E E Aft£ (*>?)( 

t,t'=±1 |fc|,|g|<7r/(2a) 



U k+q/2,r'h-q/2,r + h~q/2,T^k+q/2,T' 



X T ,T'(k,q) 



(Al5a) 
(A15b) 



E T y(k,q) 

with \ TT i(k, q) as the root of eq. flgOcft , and E TjT >(k, q) = U + r'E(k + g/2) - r£(£; - g/2), 
see eq. ( |60b| ). 

The dipole matrix elements again simplify for small q. Note that one obtains both 
contributions from q — > and g — > 7r/a. After some calculations one obtains 



|A + , + (A;;g = 0)| 2 = 



2L [3 

+ 3 



ea £7(fc) + <5?7 



(2t) 2 



ea rjE{k) — 5 



E(k)J 
(2tf 



„ „ , l2 1 / e(fc)A(jfe)(l -o" 2 



2 , 2 



(Al6a) 



(A16b) 



Note that the dipole matrix elements |A + ._(A;;g = 0)| 2 contain the contributions from q = 
it/ a for 5 = rj = 0. Eqs. (|A16|) have to be compared to the corresponding expressions for the 
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Peierls chain. It is seen that the expressions display some similarities but they show subtle 
differences. Even for q = the expressions ( |A16|) could not have been guessed. 
The corresponding dipole matrix elements become 

The dipole matrix elements between the same Peierls subbands are always strong, irre- 
spective of k or 5. However, the dipole matrix elements for transitions between different 
subbands are small for strong dimerization. Furthermore, they are small in the vicinity of 
the center and the edge of the reduced Brillouin zone. 



0)| 2 



1 

(U-2E(k)) 2 
1 

{U + 2E(k)) 2 



0)t 



(A17a) 
(A17b) 

(A17c) 



5. Spin average in the Harris-Lange model 

We need to calculate 

1 



(z^?'} = 77l ( A18 ) 
A |o> 

where 



x 



S' lr ..S' L _ 2 ^ l,r 



(5, V )x q ,(5, V )= £ 7^ E e^-^ a (l + ^6 + (-l) l (5 + V ))(l + V 5+ (~l) r (6 + V )) 

(0\S[, . . . (tzij+i _ l/tm) , S'i, ■ ■ ■ S' L _ 2 ) 



■ ■ ■ S r , (Ir+ltr Tr+llr) ? S r _ ll . . . S^O) . 

(A19) 

Since the set of spin states |0) is complete we may exactly trace it out and are left with the 
calculation of the spin matrix elements 
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M{l,r) 



- E 

2 L ^ 

S r i (Ir+ltr — Tr+llr) > S'r-D • 



(A20) 



• S[\S[, . . . S^, (hii+i - iih+i) > S'i, ■ ■ ■ S' L _ 2 ) 



The value of all spins between the sites I and r is fixed by the singlet operators at (1,1 + 1) 
and (r, r + 1). We find 



M{l,r) = 2(-l) r -'2- |r -"- a 



(A21) 



which shows that the correlation function for finding two singlet pairs at distance n — \r — l\ 
exponentially decays with correlation length £s — 1 / ln(2) . 

Since M(l,r) only depends on the distance between the two sites we may carry out one 
of the lattice sums in equation ( |A19| ). This gives 



1 



L/2-1 



n=—L/2 

n=-L/2 

The sum over n is readily taken and gives the final result 



(A22) 



71" . 



a 



7T , 



#(<?)+#(<?+-) 

a . 



(A23a) 



with the help function 



5 + 4 cos(ga) 

In particular, for the translational invariant case we find 

9(0) 



(A23b) 



2L 



(A24) 



This shows that be replaced by J g(q) / (2L) in the translational invariant case. 
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APPENDIX B: EQUATION OF MOTION TECHNIQUE FOR THE EXTENDED 

DIMERIZED HARRIS-LANGE MODEL 



Here we briefly outline the calculations for the most general case V ^ 0, 5 ^ 0. We use 
the diagonalized band picture Hamiltonian in the form of eq. (PD, and the band picture 
current operator in the form of eq. (|5%p . This has the advantage that the equation of motion 
can directly be inverted and the contribution for V = can immediately be separated. 

The equations of motions give for the four currents 

. „ , (A(u)/c)\ty + 2V(cos(ka)X° )Tl + S m(ka)X^ 
j T , T >{k, q; u) = * '- (Bl) 

where E T y = E T y(k, q) = U + r'E(k + q/2) — rEik — q/2), A r T ' = X TtT i(k, q), and X°' S T , = 
X°p(k,q). The difficult terms X c ' s = X^° + = -X^_, Y c > s = X%% = X*' + , come from the 
nearest-neighbor interaction which mixes excited pairs in the different Peierls subbands. 
With the help of eq. ([BlD, A +i+ = — A__, and A +i _ = A_ i+ we can immediately write 

(j.>o)(V) = (j UJ>0 )(V = 0) 




To determine X c,s and F c,s we have to evaluate Vu^ +q , 2 T ^k~ q /2 T \®)- Since V is simple in 
terms of the original operators Uk, h we first have to apply the inverse transformation of 
eq. (|57D . As next step we let V act and then re-transform into the operators for the Peierls 
subbands in the last step. 

The calculation shows that only four combinations of currents occur in X c ' s and Y c,s , 
namely, 

X c > s (k,q) = -f* +>+ (k,q)Jr(q) + £,_(*> q)J? e (q) (B3a) 
Y c > s (k,q) = fl_(k,q)J c { s (q) + f^ + (k,q)jr(q) (B3b) 
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with 



J i s ( ( i) = l E 



cos (pa) 



cos (pa) 



^ 2 C %) = 7 E f 

^ |p|<7r/(2a) V 



(B4a) 



(B4b) 



Since the interaction is restricted to nearest-neighbors the global, i. e., only g-dependent 
currents J\ ' 2 (q) appear in the problem. 
Equation ( p3|) becomes 

2Vf/± 



(i w >o>(v)-a>o>(v=o) = 

with GiK?) given by 

G c { s 2 (q) = - y ( cos ^ a ) 



E (?) + J?(q)G? (?)] (B5) 



|g|<7r/(2a) 



+ 



-/+.4 
/+,- 



/ 1 



1 1 

+ 



1 



uj — E-_ + u — E+ _ 



(B6) 



The quantities Gi 2 {q) are still operator valued objects since they contain \ T y. Nevertheless 
these quantities are known. We insert eq. (|59| ) and use the fact that |/ +)+ (— k, q)\ 2 = 
\f + Ak,q)\\ \UA-k,q)\ 2 = \f + ,-(k,q)\\ ™d f* +>+ (-k,q)f +> _(-k,q) = f* j+ (k,q)f +> _(k } q) 
to show that G\{q) = G c 2 (q) = 0. We set G\{q) = G x {q), G s 2 (q) = G 2 (q), = Ji(q), 

Jf(g) = J 2 {q)- They can be expressed in terms of the functions Fi 2 ^{q) of eq. ( [75[ ) as 



(B7a) 
(B7b) 



Gi(g) = itea - x+ +n/a F 3 (q) 

G 2 {q) = itea [-x+F 3 (q) + x+ +7v/a F 2 (q) 

It remains to determine J\^{q)- They can be obtained from their definitions in eq. ( p4| 
and the result from the equations of motion, eq. flBl"l), 



(B8a) 
(B8b) 



Ji(q) + ^Gt(q) = V (-J 1 (q)F 1 (q) + J 2 (q)F 3 (q)) 
Mq) + — Gt(q) = V (-Mq)F 2 (q) + J 1 (q)F 3 (q)) . 
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It is not difficult to invert these equations to obtain the currents explicitly. The result for 
the real part of the optical conductivity becomes 

Re{a(u > 0, V,5,rj)} = Re{a(u > 0,5,7/)} + (B9) 

QUJ 

1 

T\ 4- T/P, VI 4. 

-|g|<7r/(2a) 



Im i , , fl + ^(1 + FF 2 ) - (VF 3 ) 



+ 



+ V (G x G+F 2 + G 2 G+ Fx + (G x G+ + G 2 Gf)F s 



As a last step we have to factorize the products over the functions G\ y2 (q). For the 
Harris-Lange model we use eq. ( |54a|) . With the help of eq. ( |B7| ) it is not difficult to derive 
the final result for the average optical conductivity in the Harris-Lange model, eq. (173). 
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FIGURES 

FIG. 1. Band structure interpretation of the exact eigenenergies of the Harris-Lange model for 
U=2W. 

FIG. 2. Band structure interpretation of the exact eigenenergies of the dimerized Harris-Lange 
model for U=2W, 5 = 0.2. 

FIG. 3. Reduced average optical conductivity, a re d(uj > 0), in the Harris-Lange model for 
U = 2W. A broadening of 7 = 0.01IF has been included. 

FIG. 4. Reduced average optical conductivity, cr re d(w > 0,5,77), i n the dimerized Harris-Lange 
model for U = 2W, 5 = 0.2 (5 = 0.6), and rj = -0.06. A broadening of 7 = 0.01IF has been 
included. 

FIG. 5. Reduced average optical conductivity, cr rc d(w > 0, V"), in the extended Harris-Lange 
model for U = 2W and V = 0, W/2, W. A broadening of 7 = 0.01VF has been included. 

FIG. 6. Reduced average optical conductivity, a re d(u> > 0,V,6,rj), in the extended dimerized 
Harris-Lange model for U = 2W, 5 = 0.2, 77 = -0.06, V = 0, W/2, W. A broadening of 7 = 0.01W 
has been included. 
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